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Using a simple model in the context of the Dyson-Schwinger-Bethe-Salpeter approach, we investigate 
the effects of a dressed-quark-gluon vertex on psendoscalar meson masses. In particular, we focus on 
the unequal-mass case and investigate heavy-light meson masses; in addition, we study the premise 
of the effective treatment of heavy qnarks in our approach. 
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I. INTRODUCTION 

The realm of heavy-light mesons is topical and promis¬ 
ing for a number of reasons. In the light of quantum 
chromodynamics (QCD), one faces the challenge of a 
multi-scale problem, where many successful theoretical 
Ansdtze can be tested, compared, and brought together. 
With the chiral limit on one and the heavy-quark regime 
on the other end of the quark mass range, an ideal ap¬ 
proach should respect the symmetries apparent in QCD at 
both of these ends together with the other basic properties 
of the theory, e. g., its perturbative limit, and a realization 
of confinement. Overall, a nonperturbative technique is 
necessary and has the benefits of wide applicability, if 
used properly and with a model tailored to this task. 

In particular, it has been shown in [THS! that the spec¬ 
tral difference of parity partners in the sector of mesons 
with a heavy and a light valence quark is driven by a 
subtle balance between a suppressed interaction with the 
QCD ground state and an enhanced impact of light-quark 
chiral-symmetry breaking effects by the heavy-quark mass. 
The mass splitting of these parity partners is comparable 
to the ones in the sector of mesons with light valence 
quarks only, making them a suitable object for related 
investigations of dynamical chiral-symmetry breaking. 

The Dyson-Schwinger-Bethe-Salpeter-equation (DS- 
BSE) approach, a modern nonperturbative tool for quan¬ 
tum field theory imni complementary to the well-known 
lattice-regularized approach mHni, is an excellent can¬ 
didate for such a study, since all requirements are met. 
In this work we build on an earlier investigation of a 
systematic approach to dressing the quark-gluon ver¬ 
tex (QGV) and thus, consistently the quark Dyson- 
Schwinger-equation (DSE) and meson Bethe-Salpeter- 
equation (BSE) integration-equation kernels [151 El) 
which are the necessary prerequisites for a meson study. 
Vertex dressing affects the BSE kernel in the various 
mesonic channels differently, which is in accord with 
expectations from the quark model, where different terms 
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in the potential also show varying importance in different 
channels, see, e. g., El. 

In modern DSBSE studies with sophisticated effective 
interactions, a simple truncation offers the possibilities of 
comprehensive investigations, see 1181130] and references 
therein. Beyond the most popular rainbow-ladder (RL) 
truncation, there have also been sophisticated studies, 
following systematic schemata based on the structure of 
the DSE tower. However, the numerical effort [8TI - I35] 
quickly increases |36H46j and so there have been a number 
of investigations, like our present one, which make use of 
a very simple effective interaction [47] in order to be able 
to highlight particular schemes or effects, if one is able 
to sum up certain classes of diagrams or study a certain 
scheme in as great detail as possible P El HSHSi] . Apart 
from elucidating effects of high-order dressing terms or 
fully summed subclasses of diagrams, such studies provide 
a guide as to how far a systematic truncation scheme with 
a sophisticated effective interaction would have to go 
in order to achieve reliable results for various hadron 
properties. 

Our focus on heavy-light meson systems can help to 
study the importance of QGV correction terms by testing 
those in a multiply unbalanced system of one light (and at 
the same time very strongly-dressed) antiquark together 
with one heavy (and only mildly-dressed) quark in a 
relativistic bound state. Note that investigations of heavy- 
quark propagators have shown that dressing can have 
sizeable effects even for the case of b-quarks |55]-[57] . In 
contrast to earlier DSBSE studies making use of particular 
assumptions about heavy-quark propagators [5HH55] we 
would like to motivate or justify these assumptions rather 
than start out from them. Ultimately, we aim at a direct 
check of heavy-quark symmetry predictions [66] like the 
ones performed recently, e.g., in relativistic Hamiltonian 
dynamics m- 

The article is organized as follows: In Sec.|TI]we present 
the setup used for the quark DSE, the QGV, and the 
meson BSE. Results are collected and discussed in Sec. jm] 
for both the quark propagator dressing functions and the 
meson masses. After the conclusions in Sec. IIVI we have 
collected more technical details and explanations as well 
as a complete set of results in terms of tables and figures 
together with an analysis of heavy-quark effective quark 
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propagators and their comparison to the ones computed 
here in several appendices. 

II. SETUP 

Since this work is an extension of [TSl HI] , we focus on 
essentials and differences instead of repeating every detail 
here. Our calculations are performed in Euclidean space. 

A. Quark DSE 

For the purpose of investigating hadrons as bound states 
of (anti)quarks interacting via gluons, the DSE for the 
quark propagator is an adequate starting point, because 
it contains two main ingredients of such a study, namely 
the renormalized dressed gluon propagator and the 
renormalized dressed QGV r“ with the color index “. 
Its solution is the renomalized dressed quark propagator 
S for a particular quark flavor, which has the general 
structure 

=i-fpA{p^) + B{p^) (1) 

= A{p'^) y ■ p +M^)) . (2) 

In this and the following, we omit dependences on 
renormalization and regularization parameters and details 
due to the particular ultraviolet-finite model interaction 
we are using herein: the corresponding renormalization 
constants are = I. 

The dressing functions A and B characterize S com¬ 
pletely. Alternatively, one can use A and M, all of which 
are functions of the quark momentum p squared, and 
implicitly also of the current-quark mass rUg and thus the 
quark’s flavor. We write s := p^ and choose to investigate 
A(s) and M{s) below. 

The quark DSE reads 

S~^{p) = + 'B{p) = i'yp + mq + 'B{p) ( 3 ) 

= jg^Dyp-q)^^y{q)Vl{q-p) (4) 

j4 

with fg ■= f S ths quark self energy, p the quark 

momentum, and g the strong coupling constant. The 
subscript o denotes a bare quantity. With the specification 
of the effective interaction with constant strength Q to be 
used herein m via 

5^ Df,^{k) := ^(5^^ - 0^ S^{k), (5) 

one arrives at an algebraic equation for S; the same is also 
true for the other equations of relevance herein, namely 
the DSE for the QGV and the meson BSE. More precisely, 
with the explicit color prescription r“(p) = ^r^(p) and 
setting Q = 1, thereby obtaining all dimensioned quanti¬ 
ties in appropriate units of Q one obtains 

S~^{p) = i-fP + mq+-ffj,S{py{p), ( 6 ) 


where the model parameter C introduced at this point is 
defined below. At this point the quark DSE can be solved 
for any explicit form of the QGV, and results in physical 
units are obtained by providing the necessary input, i. e., 
a value for the model strength parameter Q and a quark 
mass niq for any given quark flavor. 

To proceed, one can make use of the DSE for the QGV 
in order to define a model for T^{p). In particular, again 
following [TSj, the QGV DSE is used to obtain the effective 
equation 

b) = 7 m - c 7p -S' {p) T^ b) s b) ip, (7) 

whose dependence on C stems from the effective combina¬ 
tion of the abelian and non-abelian correction terms in 
the QGV DSE such that: i) any differences in strength or 
momentum dependence of these two terms is taken to be 
proportional to the simple and herein accessible structure 
of the abelian term; ii) the sign is that of the non-abelian 
correction due to its dominant role; and iii) C is chosen 
in accordance with either lattice QCD, phenomenology, 
or another reasonable point of comparison. 

From a general point of view, the following values are of 
interest: C = —1/8 corresponds to the case with abelian- 
only dressing as, e. g., used in [49]. C — 0 corresponds to 
the popular RL truncation, which is also the zeroth order 
in the systematic schemes of the kind considered here. 
Finally, C = 0.51 was used in m as a result from fitting 
to lattice quark propagators, with good phenomenological 
success. We fix C = 0.51 throughout herein in order to 
avoid the overwhelming amounts of data that would result 
from a thorough investigation of variation in all available 
parameters. For our purposes it is best to continue from 
m for clarity of the resulting effects, easy comparison, 
and direct applicability. 

Equation 0 is the immediate basis for an iterative 
prescription, where the bare QGV serves as a starting 
value: T^ g b) = 7 m • The recursion relation is 

T^pb) = -C7p5b)r^,,_ib)^b)7p (8) 

so that at a given order n in this scheme one has for the 
QGV 

n 

= ( 9 ) 

1=0 

and the final result for the QGV is obtained by n —)■ oo. 

At this point it is also important to note that, like the 
solution of the quark DSE, also the corresponding result 
for the QGV implicitly depends on the flavor (and niq) 
of the quark it is associated with, i. e., the details and 
properties of the factors S{p) in Eqs. 0 and (]8[) . Note 
also that, since the quark-gluon interaction is not flavor 
changing, both factors S{p) in each of these equations 
must correspond to the same flavor and value of iriq. 
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B. Meson BSE 

The meson BSE reads 

[riP;k)]EF= [[K{P-,k-q)]%^ [S'(g+) r(P; q) S{q-)]cH 

Jq 

( 10 ) 

where the four-vector arguments P and k, q are the 
quark-antiquark pair’s total and relative momenta, respec¬ 
tively. E, F, G, H stand for color, flavor and spinor indices, 
T{P;q) is the meson’s Bethe-Salpeter amplitude (BSA), 
and K the fully-amputated dressed-quark-antiquark scat¬ 
tering kernel. The combination of quark propagators and 
the BSA is often combined to the so-called Bethe-Salpeter 
wave function y and one writes 

x{P;q):=S{q+)T{P;q)S{q.) , (11) 

where the meson flavor is determined by the combination 
of quark flavors from the two factors of S. The quark 
and antiquark momenta are defined as q+ = q + rjP and 
g_ = g — (1 — r])P. rj is called the momentum partitioning 
parameter, and its arbitrariness in any computation is 
equivalent to the freedom in the definition of the quark- 
antiquark relative momentum. In a covariant setup, any 
observable should be independent of ry; still, numerical 
approximations, truncations, or other model defects can 
destroy this independence, as is the case here, where 
the BSA is incomplete due to the particular form of the 
model interaction in Eq. ([^. As a result, any study using 
this particular interaction should investigate also the r]- 
dependence of numerical results; such an analysis was 
also performed in Ref. m where this particular form 
of effective interaction was introduced in RL truncation. 
Still, such a model artifact does not destroy the model’s 
capacity to elucidate our targeted meson properties and 
can be easily quantified; thus it is well under control. 

Provided one has solutions of the quark DSE for a given 
form of the quark self-energy E, taking into account a 
particular form P;^ of the QGV, one can construct a BSE 
interaction kernel consistent with the quark propagator 
dressing functions via the axial-vector Ward-Takahashi 
identity (AVWTI) [551ISS] ■ It was shown in [T^ |351135| 
that the BSE kernel corresponding to our particular setup 
when omitting contributions from gluon-unquenching [361 
155] for equal-mass constituents in the BSE reads 

T^\P-k) = -'^ J D^,{k-q)[x^X^{P-q)T,{q_,k_) 

+ 7^5(9+) Af(P;fc;q)] . (12) 

The color trace has been carried out at this point, since 
the mesonic case is straight-forward in this regard; in 
the following, we always give results with evaluated color 
traces. The superscript label ^ denotes applicability to 
a particular meson, e. g., a pseudoscalar, which we inves¬ 
tigate herein. This is important, because the structure of 
the correction term crucially depends on the structure 
of the corresponding BSA. We will detail this point below 
for the pseudoscalar case. 


While the first term on the r.h.s. of Eq. (12) is con¬ 


structed immediately from the given QGV, the construc¬ 
tion of the second term is again based on a recursion 
relation analogous to the one for the QGV. One sums 
up the correction terms up to a particular order n to get 
A^: 


Af(P;fc;9)=^A()f,(P;fc;9). 


(13) 


i=0 


The full result is then obtained by n —>■ oo. 

Generalizing to the unequal-mass case, we symmetrize 
Eq. (12) to get 


rM(p;fc)^_14y D^,{k-q)[x^x^{P-,q)V,{q_,k_) 

+ l^iS{q+) Af (P; fc; q) -f T^{q+,k+)x^[P] q)xv 
P A^J{P■k■q)S{q.)x^] . (14) 


The flavor content of this equation is encoded via the 
quark propagators and, more precisely, their arguments, 
the quark momenta k±,q±: any factor of S or with (an) 
argument (s) with subscript + corresponds to the quark 
1 with flavor 1 and those with (an) argument(s) with 
subscript _ corresponds to the antiquark 2 with flavor 
2. We note here also that in our setup for unequal quark 
masses, the heavier one is associated with the subscript 


+ ■ 

The situation becomes clearer, when we employ the 
effective interaction from Eq. ([^ to get the algebraic form 

r^(P;fc) = -^[7^X^(P;fc)rC(fc_) 

+ x^.S{k+) Af (P; k) + rC(fc+)x^(P; k)x. 

+ A^{P-k)S{k.)xu\ , (15) 

where we have written A^{P-,k) := A^{P;k;k) and 
denote the dependence on the parameter C explicitly. 
One thus finds [T5l US] : 

iA()"„(P;fc) = -7pX^(P;fc)rC_,(fc_)5(fc_)7p 

-7p^(A:+)rC_,(A:+)x^(P;fc)7p 

-XpSik+)A^^_,{P;k)Sik.hp, (16) 

where one has to take care of attributing the correct 
quark flavors and factors of S and Ej, according to their 
arguments with respect to the subscripts ±, as described 
above. This expression can be computed, the recursive 
terms for A^^{P;k) summed up to any desired order, 
and the resulting algebraic equations can be solved. It is 
helpful to note here that also for n —>■ oo the equations 
can be solved explicitly via geometric summation. 

The initial condition for this recursion relation is [49] 


A“(P;fc)=0, 


(17) 


which (only) in the pseudoscalar case for equal-mass 
quarks and 7 = 1/2 implies [15] 

A^^(P;fc)=0 ^ A^’^(P;fc) = 0, (18) 
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FIG. 1. Dressing functions A(s) and M{s) for different values of n, corresponding to a quark mass m„ = O.OlGeV. 


which is an exceiient testing case for our generai setup. 
Further detaiis on this construction are rather technicai 
and are thus coiiected in App. 

III. RESULTS AND DISCUSSION 

In the present investigation we ask the question how 
important corrections to the popuiar RL truncation are 
in the case of mesons with unequai-mass constituents. 
Whiie our resuits aiiow not oniy quaiitative, but aiso 
quantitative statements, caveats are due to the simpiic- 
ity of the interaction and the resuiting oversimpiification 
of structures reiated to reiative quark-antiquark momen¬ 
tum, as weii as the resuiting artihciai dependence on the 
momentum-partitioning parameter 77 . As a resuit, we not 
oniy attempt to study and quantify effects reiated to the 
number of recursion steps used for the QGV and the BSE 
kernei, which is our main objective, but aiso investigate 
the dependence of reievant quantities on 77 . The subsec¬ 
tion on the quark DSE resuits, however, is independent 
of 77 . 

Our modei parameters are fixed to: C = 0.51 sets the 
amount and quaiity of dressing in the DSE of the QGV, 
and is fixed throughout this paper, as mentioned above. 
Investigations of the effect of using different values of C 
have been performed in [15] for equal-mass constituents to 
some extent, but mainly concerning the change from the 
abelian value of C = —1/8 to a non-abelian-dominated 
version, which is phenomenologically reasonable and ex¬ 
tendable. A detailed investigation of the C-dependence of 
our results is beyond the present study and thus deferred 
to subsequent publications. 

The effective coupling is set toG = 0.69, a value compat¬ 
ible and preadjusted to fit vector-meson masses through¬ 
out entire quark-mass range. This fact also explains why 
our results for pseudoscalar mesons as they are presented 
below in Fig. do not appear completely satisfactorily on 
the level of a pure theory-experiment comparison. This 
is, however, not the main point of our study and serves 
mainly to keep all masses and quantities under investiga¬ 


tion on their respective domains of reasonable values. 

Finally, we also take the current-quark masses directly 
from [15] and set niu = O.OlGeV, ms = 0.166GeV, me = 
l.SSGeV, mi, = 4.62GeV. For the pion we assume isospin 
symmetry and the equality of the current-quark masses 
of the u and d quarks. 


A. Quark DSE 

The equations resulting from the general setup de¬ 
scribed in Sec. [ill after applying the effective interaction 
defined in Eq. , are algebraic and can be tackled via 
standard root-finding techniques. We note here that the 
solutions from equations involving polynomials can be 
defined piecewise, as is also the case in our investigation 
and visible in Eigs.[^and|^ 

A typical set of solutions for the dressing functions A(s) 
and M{s) is plotted in Fig. [^for the light quark mass 
and different steps in the recursion, starting from n = 0, 
which corresponds to RL truncation, and including the 
fully dressed vertex represented by 77 = 00 . Analogous 
plots for the strange, charm, and bottom masses can be 
found in Fig. in App. [Aj 

It is notable, how the simple behaviour of A and M 
are changed from the RL result (dashed lines) with the 
QGV dressing switched on. The change is more drastic 


TABLE I. Galculated mass function M{s) at s 
and heavy quarks, in GeV. 

= 0 for light 


M„(s = 0 ) 

o' 

II 

Meis = 0 ) 

Mt{s = 0 ) 

n = 0 

1.011 

1.182 

2.556 

6.972 

n = 1 

0.731 

0.941 

2.481 

6.966 

n = 2 

0.837 

1.027 

2.497 

6.966 

n = 3 

0.798 

0.998 

2.494 

6.966 

n = 4: 

0.814 

1.009 

2.494 

6.966 

n = 00 

0.810 

1.006 

2.494 

6.966 
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FIG. 2. n-dependence of various quantities: dressing functions M{s) and meson masses, all in GeV. Left panel: M(s = 0) 
for all four quark masses; center panel: Calculated mass function M{s) at s = —For the unequal-mass case the mass 
function of the heavier quark is presented. The experimental values are used for the meson masses mn (see text). The labels H 
of the subplots denote n, K, D, B; right panel: Pseudoscalar bound-state masses for the equal-mass case calculated from the 
BSE for different quark flavors. 


for A, where one can find qualitatively different trends 
on the timelike domain at the scales relevant for the 
respective bound states, i.e., s = — m^/4, where ttim 
is the bound state’s mass. For M the changes are less 
pronounced; however they are certainly quantitatively 
relevant. In addition, M appears together with A in 
the denominator of the dressed propagators, where their 
effects are combined. 

To illustrate this further, we present both tabulated 
concrete values as well as plots. First of all, we tabulate 
the mass function M(s) at s = 0, which is one possible 
definition of a constituent-quark mass, in Tab. |I] The 
same information is plotted in the left panel of Fig. 
One observes two main features: First, the results seem 
to alternate in the sign of the effect relative to the fully 



FIG. 3. Relative differences Rn of dressing functions M(0) 
for n loops compared to the case n —>■ oo, for all four quark 


dressed result with respect to n: odd n yields a lower 
result than n —>■ oo, even n a higher one. Secondly, as 
expected, dressing effects are weaker in the heavy-quark 
domain. While the changes are pronounced in the light- 
quark case, a result identical to the one for n —> oo is 
reached already for n = 3 for the charm quark, and n = 1 
for the bottom quark, respectively. Defining the relative 
changes from RL truncation {n = 0) to the fully dressed 
result (n —)• oo) as 

M„(s = 0)-Moo(s = 0) 

■= -M„(, = 0)- ■ <“> 

we obtain i?o-values of 0.25 for the u quark, 0.17 for the s 
quark, 0.025 for the c quark, and 0.00086 for the b quark. 
This comparison, together with values of for all n and 
quark masses are plotted in Fig. 

In a similar fashion, we present values for the mass 
functions for different quarks at that value of s which 
appears as the argument of the dressing functions of S 
in the solution of the BSE. For the case of the effective 
interaction ([^ used herein this is proportional to rf' and 


TABLE 11. Calculated mass function M(s) at s = 
in GeV. For the unequal-mass case the mass function of the 
heavier quark is presented. The experimental values are used 
for the meson masses mu (see text). 


H 

TT 

K 

D 

B 

be 

Vb 

n — Q 

1.014 

1.204 

2.620 

7.013 

2.733 

7.157 

n = 1 

0.731 

0.943 

2.516 

7.004 

2.559 

7.123 

n = 2 

0.839 

1.041 

2.540 

7.004 

2.613 

7.126 

n = 3 

0.799 

1.005 

2.535 

7.004 

2.598 

7.126 

n = 4: 

0.815 

1.019 

2.536 

7.004 

2.603 

7.126 

n = OO 

0.811 

1.015 

2.536 

7.004 

2.602 

7.126 


masses. 
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the pseudoscalar bound-state mass-squared. While we 
investigate the entire range of values of ry from 0 to 1 
below, we use rj = 0.5 for the purpose of this argument 
and set s = —171^/4. The corresponding numbers are 
collected in Tab. mi and the first four columns are also 
plotted in the center panel of Fig. For the unequal- 
quark-mass case we give/plot values of the mass function 
corresponding to the heavier quark flavor. In order to 
avoid additional uncertainty, the experimental values are 
used for the meson masses niH for each meson as given in 
the column labels: = 0.140 GeV, mx = 0.494 GeV, 

niD = 1.870 GeV, tub = 5.279 GeV, = 2.984 GeV, 

and rurj^ = 9.398 GeV. It should be noted here that in 
the heavy-light case with values of 77 ranging from 0 to 
1 one potentially covers the whole range s = —m'j^ to 
s = 0, corresponding to ?y = 1 and 77 = 0, respectively. 
This is particularly interesting regarding the behavior of 
the dressing functions with n as discussed in more detail 
in App. 1^ 

The pattern in Tabs. |T] and |TT] are similar, which confirms 
a consistent picture of the dressing effects on the functions 
A(s) and M(s) on the entire relevant parts of the timelike 
s domain. This remark remains also valid if one uses a 
more sophisticated effective interaction than the one of 
Eq. (§: in such a case, the integral-equation character 
of the BSE remains intact and the domain on which A(s) 
and M{s) are sampled most prominently in the complex 
s plane is the one surrounding the negative real axis, 
reaching out to a scale of s = —rri\j/4\ for more details, 
see the appendices of [22l[27l[70H72]. 

It is illustrative to note here that the alternating conver¬ 
gence pattern between even and odd n can be attributed 
to the negative sign of the r.h.s. of Eq. ([^, the recursion 
relation for the QGV. 


B. BSE 

With the details of the solutions of the quark DSE 
laid out and the solutions obtained for the four relevant 
quark flavors, we proceed to the corresponding results of 
the pseudoscalar meson BSE. As a first test and for easy 
comparison and anchoring, we set 77 = 1/2 and recompute 
the equal-quark-mass results of Ref. m, which we have 


TABLE III. Pseudoscalar bound-state masses for equal quark 
masses calculated from the BSE for different quark masses, in 
GeV. 



TT 

07 + 


Vb 

n = 0 

0.149 

0.669 

3.033 

9.505 

n = 1 

0.132 

0.604 



n = 2 

0.140 

0.639 

3.012 

9.504 

n = 3 

0.137 

0.626 



n = 4: 

0.138 

0.632 

2.998 

9.500 

n = 00 

0.138 

0.631 

2.985 

9.489 


TABLE IV. Bound-state mass for the kaon as a function of n 
and ? 7 , given in GeV. 


V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

n = 0 

0.455 

0.455 

0.458 

0.461 

0.464 

0.472 

0.483 

n = 1 



0.399 

0.409 

0.416 

0.417 


n = 2 

0.421 

0.426 

0.433 

0.436 

0.440 

0.447 

0.453 

n = 3 


0.399 

0.420 

0.426 

0.431 

0.431 

0.427 

77 = 4 

0.407 

0.416 

0.426 

0.431 

0.435 

0.441 

0.443 

n = 00 

0.392 

0.411 

0.425 

0.430 

0.434 

0.439 

0.438 


collected in Tab. and plotted in the right panel of Fig. 
While the emerging convergence pattern is qualitatively 
similar to the ones observed above for the quark mass 
function, we also encounter cases where the BSE does not 
yield a solution for particular values of n, e.g., n = 1,3 
for the 77 c and rjb masses. 

To understand this, a few remarks are in order. Eirst 
of all, the homogeneous BSE is an equation, which in 
principle does not need to have a solution in all possible 
cases or for all given configurations in a certain situation or 
setup. Depending on truncation setup, model features or 
defects and also ranges of parameters, one can encounter 
cases where there is no solution. In such a case, we denote 
the missing solution by three dots in the corresponding 
spot in the table and curves in figures go across such a 
point. In Fig. for missing results in a particular case, 
no calculated data point is generated. 

That said, we proceed to the case of unequal (anti)quark 
masses in the meson and present results for the various 
combinations of quark flavors. We start with the kaon, for 
which the results are presented in Tab. |IV| At this point 
a comment about the 77 dependence of the results is neces¬ 
sary. Eq. leads to the model artifact that the relative 
momentum in all parts of the calculation of the meson 



FIG. 4. Bound-state mass for the kaon as a function of n and 
rj, given in GeV. Even n are depicted by dashed lines, odd 
ones by dotted lines, and the fully summed result by a solid 
line. 
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mass vanishes. This not just simplihes the structure of 
the equations but also the structures of the pseudoscalar 
(and other) BSAs, the QGV, and the correction terms A. 
One of the results of these simplifications is the artificial 
77 dependence of the results already mentioned above. To 
remain in control of our results and the corresponding 
discussion, we thus have to analyze and quantify the 77 
dependence of the observables under consideration. We 
do this by providing a complete set of results as well as 
encoding the 77 dependence in the form of a systematic 
error in our results in Fig. the error bars are plotted 
from the lowest to the largest value of the mass result for 
a given n. As a result, they are asymmetric and the value 
of 77 chosen for the data point, as defined below, may on 
occasion be also either the smallest or the largest value 
available at this n. 


To illustrate the extent of this effect, we present a con¬ 
crete example, namely the kaon mass as a function of both 
n and 77 tabulated in Tab. |IV| As a further illustration, we 
present the same numbers plotted in Fig.[^ One observes: 
convergence with n; an alternating convergence pattern 
with respect to odd and even n analogous to the ones 
observed above; and a moderate dependence on 77 in the 
sense that the effect across the entire rj range is compara¬ 
ble in relative size to the effect of the dressing from RL 
truncation to the fully summed result for a fixed value of 
77 . Concretely, we observe changes of the order of 10%, 
which is in agreement with the analysis in 03, where the 
authors quote changes smaller than 15%. Judging from 
the particular behavior of the curves in Fig. at the outer 
boundary of the 77 G [ 0 , 1 ] interval, one can see that these 
extreme values of 77 also seem to lead to correspondingly 
extreme values of rriK- While it is certainly correct to 
state that the error bars in Fig. should represent the 
entire range of 77 , it is also fair to remark that in practice 
the extreme 77 values may not be representative to an 
amount that actually justifies the size of these error bars 
and we in general regard them as overestimates of more 
suitably defined systematic errors. 


Tables and figures analogous to Tab. W and Fig. 
are collected—for all mesons considered here and whose 
masses are presented below in Fig. [ 5 ]— in App.[^in Tab. |VI| 
and Fig.[^ There, we have several rows completely filled 
with dots indicating that no solutions exist for a given n. 
However, there may still be small domains in between the 
values of 77 listed in the figures. While we did not include 
those in Tab. VI and Fig.[^to remain at a comprehensible 
set of results, we have done a finer search and included 
them in the results presented in Fig. These cases are 
easily recognized by their small error bars, which we chose 
not to rescale or blow up artificially; a prominent example 
is the n —>■ 00 value for the r\c mass in the lower right 
corner of the figure. 


Finally, in Fig. we collect all our results in a com¬ 
pact fashion and compare them to experimental data. 
As already mentioned above, the aspect of comparison 
to experiment is not central in our argumentation and 
serves merely to put our results into perspective. What 


TABLE V. Absolute and relative mass differences for the pion 
and the unequal-mass pseudoscalar mesons with all possible fla¬ 
vor combinations. Amu is given in GeV, the other quantities 
are dimensionless (see text). 


H 

AniH 

Amjj' 

TV 

0.011 

0.078 

Vc 

0.048 

0.016 

lib 

0.016 

0.002 

K 

0.031 

0.072 

D 

0.059 

0.034 

B 

0.124 

0.025 

Ds 

0.074 

0.039 

Bs 

0.099 

0.019 

Be 

0.150 

0.024 


we focus on here are the effects of the dressing introduced 
in the QGV for the unequal-mass case as compared to 
the analysis already provided in [15]. In Fig.j^we present 
pseudoscalar ground-state masses for the pion as a ref¬ 
erence and all flavored pseudoscalar mesons up to and 
including b quarks. The respective experimental values 
are marked by the horizontal lines in each of the subplots, 
while our calculated results are given as filled circles with 
error bars as discussed above. In particular, the filled 
circles are obtained via the following 77 values: 0.5 for the 
TT, K, 77 c, and r]b, 0.82 for the D, 0.91 for the B, 0.8 for 
the Ds, 0.89 for the Bg, and 0.79 for the Re- 

Clouded to some extent by the systematic errors, we 
can still observe dressing effects in any given 77 slice of our 
data. Using the values plotted as the filled circles in our 
figure, we calculate absolute and relative changes from 
the dressing effects in a comparison of the fully dressed 
result to RL truncation as follows: the absolute difference 
in mass mn for a given meson H is obtained between 
fully dressed and RL result and denoted by Amn- The 
corresponding results for our set of states depicted in 
Fig. is tabulated in the hrst data column of Tab. |V| 
and given in units of GeV. The second column lists the 
corresponding relative difference 


Am5|' := 


Ti —n 

— m 


H 


( 20 ) 


‘■H 


which is dimensionless and a better point of comparison 
to, e. g., our estimates of the systematic error from our 
model artifacts. While we compare the difference to the 
fully dressed result here, one may also divide by the RL 
result instead, which can be uniquely related to Eq. (20). 


Since the differences are small, however, the effect of such 
a change does not affect our discussion. 

The absolute difference is largest in the unequal-mass 
case involving a b quark, topped by the Be mass, which 
changes by 150 MeV from RL truncation to the fully- 
dressed case in sharp contrast to pseudoscalar bottomo- 
nium. While the corresponding relative differences of « 2 












FIG. 5. Bound-state masses for tt meson, r/c, and flavored psendoscalar gronnd states as a function of n, given in GeV. The 
dependence on rj is illustrated via the error bars. Galculated results are given by blue dots; experimental data are represented 
by horizontal lines. 


% make this effect look small again, one must not forget 
that this is still much larger than relevant spectroscopic 
properties like, e. g., the hyperfine splittings between the 
0 and 1 ground states. In order to concretely inves¬ 
tigate the dressing effects on splittings, we need a similar 
result like the one presented here for pseudo scalars in 
addition also for the vector case. While this is beyond 
the scope of the present study, it will be investigated and 
presented in future publications. 

Regarding the relative differences, we observe that for 
our data points their size is similar to the equal-mass case 
and of the same order of magnitude throughout. The rj^ 
and rib relative mass differences are at the lower end of the 
set of results presented here. In fact, this is also true if one 
considers the absolute mass differences, where the effect 
tor rjb is comparable to the tt and for rj^ is comparable to 
the K. 

It is important to note at this point that the pion 
results are somewhat influenced by the pion mass’ chiral- 
limit behavior. In particular, the AVWTI is satisfied at 
every recursive step n in our setup and thus in each case 
the pion mass is zero. Therefore, the variation of the 
pion mass with n is expected to be small from the very 
beginning. The typical size of the mass difference of « 10 
MeV confirms this expectation. 

Overall, we see that the effects are not negligible, 
warranting further systematic exploration of effects on 
beyond-RL treatment of meson (and baryon) properties 
in the DSBSE approach. At the same time, our results 
provide no reason to disregard careful and comprehensive 
studies in RL truncation from the outset, which opens 
the door for straight-forward and large-scale phenomeno¬ 


logical studies in an RL setup. 


IV. CONCLUSIONS 

We have extended a previously and phenomenologi¬ 
cally successfully applied systematic DSBSE truncation 
scheme and setup for pseudoscalar mesons from the equal 
to the unequal-mass case for the first time. Constructing 
and reviewing the necessary pieces in the pseudoscalar 
quark-antiquark BSE, we have studied and quantified 
effects arising from simplification artifacts in the model 
used. Still, justification for using such a simple model 
comes from the fact that it makes involved investigations 
of truncations and/or other features within the DSBSE 
approach feasible. Our results are presented in a com¬ 
prehensive manner in both tables and figures for easy 
comparability and straight-forward discussion. Under the 
assumption that the rj dependence apparent in our results 
can be quantified and treated as a systematic error, we 
have presented meson masses for RL truncation as well 
as for several steps in the truncation scheme, including 
the fully summed case. 

The results show clear and understandable patterns and 
confirm important and central expectations according to 
the symmetries at both the chiral and the heavy-quark 
ends of the quark-mass range available from experimental 
data. Dressing effects compared to RL truncation are 
sizeable, but not overwhelming, which provides support 
for both RL studies targeted at phenomenology as well as 
investigations of hadron properties beyond RL truncation. 
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FIG. 6. Quark propagator dressing functions A(s) and M{s) as functions of n, where n = 0 and n = oo are represented by the 
dashed and thick solid lines, respectively. Upper left and right panels: >l(s) and M{s) for nia = 0.166GeV; middle left and right 
panels: A(s) and M{s) for rric = l.SSGeV; lower left and right panels: A(s) and M{s) for mi, = 4.62GeV. 


Our tests of the validity of an effective Ansatz for the 
quark propagator using constant dressing functions adds 
to the consensus about the bottom quark being essentially 
heavy, while the charm quark is not. 
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Appendix A: Quark propagator dressing functions 

In this appendix we present figures of the dressing 
functions A(s) and M{s) for the s, c, and b quarks, anal¬ 
ogous to Fig. [^presented in the main text. This serves 
to provide a complete set of results together with the 
corresponding illustrations. All dimensioned quantities 
are given in GeV. 

One observes easily how the distinctive features of the 
curves move further into the timelike domain with in¬ 
creasing quark mass. In particular, a scale relevant to the 
computation of meson properties via the BSE in our setup 
is the value s = —mjj/A for which M{s) is also tabulated 
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FIG. 7. Comparison of j4(s) for the four quark masses, plotted on the same domain for each n. Left upper panel: n = 0; right 
upper panel: n = 1; Left middle panel: n = 2; right middle panel: n = 3; Left lower panel: n = 4; right lower panel: n = oo. 


in Tab. [TI] While for the light mesons the meson mass 
range up to 2 GeV implies values of s « — 1 and higher 
as the region of interest, the corresponding range for bot- 
tomonium would be s « —25 and higher. Apparently, 
said features in A and M as well as the regions where the 
differences among the curves resulting from the various 
values of n are most pronounced are relevant regarding the 
BSE for light quarks and become unimportant for heavier 
quarks. To elucidate this further, we provide a detailed 
comparison to effective free model-quark propagators in 
App.0 

At this point it is also interesting to note that, as already 
indicated in Sec. m depending on the value of ij one 
potentially covers the whole range of arguments s = 


(for 77 = l)tos = 0 (for ?; = 0) in the quark propagators 
as they appear in the BSE. In practice this means that one 
can attempt to choose r] such that the distinct structures 
in the respective dressing functions are avoided by the 
interplay of s+ = —rfrri^j vs. s_ = —(1 — ri)^rrijj, thus, 
e.g., minimizing dressing effects in the resulting meson 
properties. 


Appendix B: Heavy-quark propagator 

In this appendix we investigate and illustrate the prop¬ 
erties of the quark-propagator dressing functions A and 
M for heavy quarks. This allows us to provide support 












11 




^ ^ % 

^ \ 

;__\| 


--- mu jl 

ms 

me / / 

’ - - nib 

. 

* ^ 

*'•’.. 


-25 -20 -15 -10 -5 0 5 10 


A 

1 1 
/ 1 
/ \ 

y \ 

* 

1 ■ ■ - nau 

1 ---- ms 

1 - - 

3 - - m/j 




-25 -20 -15 -10 -5 0 5 10 


\ it 

'v*'' 

- 

•_--- niu 

; ---- ms 

--me 

: - - nib / 

. 

\ "".. 


-25 -20 -15 -10 -5 0 5 10 

■ !; 

'■ h 

— ms '' S 

® II S 

f » I'V 

--me / \ :r. 

/ \ mV 

-• nib / at 

- 

* M I 

^ * , 

* 

. .. 

■.....;■ 


-25 -20 -15 -10 -5 0 5 10 


FIG. 8. Comparison of the relative difference of our calculated A{s) to an Ansatz j4(s) = 1 for the four quark masses, plotted 
on the same domain for each n. In particular, we plot A(s) — 1. Left upper panel: n = 0; right upper panel: n = 1; Left middle 
panel: n = 2; right middle panel: n = 3; Left lower panel: n = 4; right lower panel: n = oo. 


for Ansdtze used for the quark propagator in the heavy- 
quark case, where A(s) is approximated by 1 and M{s) 
is approximated by a suitable value Mg, both of which 
are assumed constant. 

To check this approximation we compare our calculated 
propagator dressing functions to the assumed effective 
version with A(s) = 1 and M{s) = M{0). As a first 
illustration we present A(s) for all four quark flavors un¬ 
der consideration on the same domain, namely from the 
relevant timelike point in bottomonium to the spacelike 
region. In Fig. we plot the four curves correspond¬ 
ing to the four quark masses separately for each n and 
find qualitatively interesting features: First of all, on the 


considered domain, for all n, the behavior of A(s) resem¬ 
bles that of a constant only for the b quark; the c-quark 
propagator is still considerably dressed. Secondly, while 
one has to be aware of relevant domains for the lighter 
quarks regarding their respective bound-state-mass de¬ 
fined scales, the first point remains valid also if one only 
looks at an appropriately rescaled smaller portion of the 
timelike domain. 

As a next step, we compare both calculated dressing 
functions A and M to their effective constant counterparts 
by plotting the relative difference between dressed and 
constant versions. The results for A(s) are shown in 
Fig.|8l again each group of the four quark mass versions 
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FIG. 9. Comparison of the relative difference of our calculated M{s) to an Ansatz M{s) = const. = M(s = 0) for the four 
quark masses, plotted on the same domain for each n. In particular, we plot [M{s) — M(0)]/M(0). Left upper panel: n = 0; 
right upper panel: n = 1; Left middle panel: n = 2; right middle panel: n = 3; Left lower panel: n = 4; right lower panel: 
n = 00 . 


separately for the usual values of n. The same is shown 
in Fig. 1^ for M{s). For the constant effective value of 
Mq we choose M{s = 0), wich is one possible definition 
and for our purposes is as good as any other choice on 
that domain of s on which dynamical chiral symmetry 
breaking is clearly visible. 

Apparently, our figures suggest that the transition from 
dressed to effective quark propagator happens somewhere 
between the c and b quark masses, which is not entirely 
unexpected: charm quarks are in general not regarded 
as heavy in the sense that effective theories work there 
without trouble; bottom quarks, on the other hand, are 


usually regarded as well-approximated by an effective 
form of the quark propagator, which we find as well. The 
key in Figs. and in order to decide whether or not one 
is close to the form of a free propagator is the proximity 
of the relative difference to zero and a constant type of 
behavior on the relevant domain. The only case in our 
computed results where this is clearly visible, is the case 
of the b quark. 
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FIG. 10. Meson bound-state masses as functions of n and rj, given in GeV. Even n are depicted by dashed lines, odd ones by 
dotted lines, and the fully summed result by a solid line. If no solution is found, no surface is plotted at the corresponding n. 
Left upper panel: pion; right upper panel: D; Left upper center panel: B-, right upper center panel: Ds', Left lower center panel: 
Ba', right lower center panel: Be Left lower panel: rje', right lower panel: rjb- 
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Appendix C: jy-dependence of meson masses 


In this appendix we collect data and plots about the 
details of the dependence of the meson masses on the 
momentum-partitioning parameter 77 as a function of the 
order n in our scheme in analogy to Tab. IV and Fig. 
in the main text of the paper. The data is collected in 
Tab. |VI| the corresponding plots are shown in the various 
panels of Fig. As a reference and excellent example 
for the kind of behavior to expect, we present the pion, ijc, 
and rjb cases in addition to the unequal-mass collection 
of data and plots. For the pion the alternating pattern 
of convergence of the odd and even n numbers described 
above in connection with the kaon case is reobserved. For 
the other cases, the situation is somewhat obscured by the 
lack of solutions for odd n on our main 77 grid points. Still, 
one can see convergence with n as well as a pronounced 
77 asymmetry for the heavy-light case, which is the source 
of the large error bars plotted in Fig. Notwithstanding 
this, the results clearly corroborate the systematic char¬ 
acter of both the approach and the truncation scheme 
presented here as well as the validity of qualitative as well 
as quantitative statements made above. 


Appendix D: pseudoscalar kernel details 


Following Refs. [15] [49] we start from the recursion 
relations for the QGV F^, Eq. (| 8 | and the BSE correction 


term Eq. (16). The first step in our construction is 

the decomposition of each of these quantities in terms of 
Dirac covariants. 

The full QGV has 12 covariant structures built from 7 ^ 
and the two independent four-vectors of the quark and 
antiquark lines. After the effective interaction, Eq. 
has been employed, only three of those are nonzero and 
one obtains 


F^(p) = 01(77^)7^ -t a 2 {p^)l -pPbi - ia^{jp^)p^ 


(Dl) 


With this structure and the initial condition that the 
QGV be bare 


= 7m 


(D2) 


one can express the QGV via its recursion relation to any 
desired order in terms of the covariants given in Eq. (Dl) 


and the quark propagator S{p). The result, in turn, can be 
inserted in the quark DSE and yields algebraic equations 
for A(s) and M{s) via Dirac-trace projections onto the 
two covariant quark propagator structures. 

A similar strategy is used to compute Am^(T’). One 
starts off with finding a suitable decomposition in terms 
of Dirac covariants for the quantum numbers appropri¬ 
ate for the meson M under consideration, in our case 
pseudoscalar. In our setup the pseudoscalar BSA has 2 
nonzero components from the general four: 

Tps{P) = fiiP^) ^lB + f 2 {P^)lBl■P (D3) 
with the unit vector P := 


The corresponding Aq-,^ has in general 12 covariant 
structures, four of which are nonzero in our particular 
setup. We have, omitting the subscript denoting the 
pseudoscalar case 

A^(P) = /3i(p2) 757^ + ^2(p2) 75 7 • P Pm 

+ ^3(P")754 + /^4(P")75 7m7-P- (D4) 


One can obtain the scalar functions j3 := {/3j}, j = 
1,2, 3,4, from a recursion relation extracted from the 
recursion relations for the QGV and A via defining the 
projection operators 

Pm.i := ^7m 75 (D5) 

Pm ,2 := ^Pm7-P75 (D 6 ) 

Pm,3 := ^Pm 75 (D7) 

Pm,4 := ^ • P 7m 75 • (D8) 

Then at every i one has 

P] = (Xps),feTr [Pm,7cA;,p,] , (D9) 


where the matrix A4ps is given by 


/ 


•^^ps — 


- -- 0 

3 3 

_4 4 0 

3 3^ 

0 0 I 

V 0 0 -i 


0 \ 
0 


/ 


(DIO) 


and one can write (see also Ref. [49] for more details) 


/3* = Mp,r. 


(Dll) 


The vector P*, in turn is obtained, writing a := { 07 }, 
j = 1, 2, 3, from 


f* = G_dr^-f G+d|"^-f (D12) 

which can be calculated when the matrices G and L are 
known. d+ and a_ denote the coefficients of the QGV 
decomposition corresponding to the + and _ arguments 
appearing in their dehning quark propagators as given 
above. Similarly, we define 


B_ 

:= M(p^_)A(p^_) 

(D13) 

B+ 

:= M(pI)A(pI) 

(D14) 

A_ 

:= (77 — 1 ) Vp^A{p^) 

(D15) 

A+ 

:= 77 '/l^A{p'^) 

(D16) 

A_ : 

= B^{p^_) + P^A^{p^_) 

(D17) 

A+: 

= B^{pI)+P^A^{pI) 

(D18) 


and obtain 
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TABLE VI. Bound-state masses for various mesons the as functions of n and rj, given in GeV. 
TT D 


n 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

n = 0 

0.1490 

0.1488 

0.1486 

0.1486 

0.1486 

0.1488 

0.1490 

n — 9 

1.645 

1.651 

1.673 

1.691 

1.716 

1.797 

1.946 

n = 1 

0.1309 

0.1316 

0.1319 

0.1320 

0.1319 

0.1316 

0.1309 

n — 1 








n = 2 

0.1397 0.1398 

0.1398 

0.1398 

0.1398 

0.1398 

0.1397 

n — 2 

1.763 

1.731 

1.709 

1.702 

1.701 

1.749 

1.900 

n = 3 

0.1364 

0.1368 

0.1370 

0.1370 

0.1370 

0.1368 

0.1364 

n — 3 








n = 4 

0.1378 

0.1381 

0.1382 

0.1382 

0.1382 

0.1381 

0.1378 

n = 4 

1.795 

1.752 

1.718 

1.705 

1.696 

1.731 

1.881 

n = oo 

0.1374 

0.1377 0.1379 

0.1379 

0.1379 

0.1377 0.1374 

71 = 00 












B 







Ds 





V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

V 

0 

0.2 

9.4 

0.5 

0.6 

0.8 

1 

n = 0 

4.801 

4.807 

4.825 

4.842 

4.868 

4.983 

5.456 

n — 0 

1.815 

1.819 

1.839 

1.857 

1.881 

1.954 

2.076 

n = 1 








n — 1 








n = 2 

5.017 

5.004 

5.003 

5.006 

5.012 

5.034 

5.370 

n — 2 

1.933 

1.896 

1.870 

1.862 

1.861 

1.911 

2.047 

n = 3 








n = 3 








n = 4 

5.069 

5.050 

5.041 

5.040 

5.041 

5.046 

5.350 

n = 4 

1.960 

1.915 

1.877 

1.862 

1.852 

1.894 

2.036 

n = OO 








71 = 00 






1.880 






B 

S 







Sc 





V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

n = 0 

4.960 

4.965 

4.984 

5.002 

5.029 

5.150 

5.567 

n — 0 

6.147 

6.150 

6.172 

6.193 

6.226 

6.329 

6.471 

n = 1 








n — 1 








n = 2 

5.171 

5.156 

5.154 

5.156 

5.161 

5.181 

5.497 

n — 2 

6.325 

6.299 

6.286 

6.283 

6.281 

6.319 

6.507 

n = 3 








n — 3 








n = 4 

5.150 

5.199 

5.189 

5.187 

5.187 

5.189 

5.483 

n = 4 

6.361 

6.327 

6.306 

6.298 

6.290 

6.312 

6.521 

n = OO 








71 = OO 




















Vb 





V 

0 

0.2 

0.4 

0.5 

0.6 

0.8 

1 

V 

0 

0.2 

9.4 

0.5 

0.6 

0.8 

1 

n = 0 

3.062 

3.037 

3.033 

3.033 

3.033 

3.037 

3.062 

n — 9 

9.487 

9.479 

9.498 

9.505 

9.498 

9.479 

9.487 

n = 1 








n — 1 








n = 2 

3.161 

3.082 

3.024 

3.012 

3.024 

3.082 

3.161 

n — 2 

9.607 

9.555 

9.521 

9.504 

9.521 

9.555 

9.607 

n = 3 








n = 3 








n = 4 

3.185 

3.090 

3.018 

2.998 

3.018 

3.090 

3.185 

n = 4 

9.629 

9.566 

9.524 

9.500 

9.524 

9.566 

9.629 


n = oo ... ... ... 2.985 ... ... ... n = oo ... ... ... 9.489 


L = 


-C 

A_ A+ 


/ 2B_B+- A_A+ 
2{B_B++A_A+) 
4i (A_B+ - A+B^) 
V i(A_B+ - A+B_) 


i (B_B++A_A+) 
2(B^B++A_A+) 
4i (A_B+ - A+B_) 
i(A_B+ - A+B_) 


^i(A_B+-A+B_) 
2i (A_B+ - A+B_) 
-4 (B_B+ + A_A+) 
-(B^B++A_A+) 


i(2A_B++A+B_) \ 
2i (A_B+ - A+B_) 
-4{B_B++ A_A+) 
-{B.B++A_A+) 


(019) 


The two matrices G_ and G_|_ are associated with the corresponding quark propagators with the _|_ and _ arguments 
as defined above: 
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G_ = 


-C 

A2 A+ 


f i {B+Al + B_A+A_ + 2BlB+) /i + i {A+Al - B_B+A_ + 2BlA+) /z 
2i {-B+Al + 2B_A+A_ + B‘iB+) /i + 2i {-A+Al - 2B_B+A_ + BlA+) /z 
-4 {A+At + 2B_B+A_ - BlA+) /i - 4 {-B+A‘i + 2B_A+A_ + BlB+) /a 
V {-A+Al_ - 2B_B+A_ + BlA+) /i + {B+A'i - 2B_A+A_ - BlB+) /a 


(^B+A'i - 2B_A+A_ - BlB+) fi{r] - l)^ - {A+A'i + 2B_B+A_ - BlA+) /2(r; - l)^ 
-2iP^ {B+A'i - 2B_A+A_ - BlB+) /i(?? - l)^ - 2iP'^ {A+A^_ + 2B_B+A_ - BlA+) /2(?? - 1)^ 
-4P2 (A+Al + 2B_B+A_ - BlA+) fi{r] - l)^ - 4p2 {-B+At + 2B_A+A_ + B'iB+) /2(r? - 1)^ 
-p2 [A+Ai + 2B_B+A_ - BlA+) fi{r] - l)^ - P^ {-B+Al + 2B_A+A_ + B‘iB+) f2{r] - if 


-fVP^iv - 1) {A+Al + 2B.B+A_ - BlA+) f - \i^fP^{ij - 1) {-B+A^_ + 2B_A+A_ + BlB+) /2 \ 
-2iVP^{r] - 1) {A+Ai + 2B_B+A_ - BlA+) f - 2iVP^{ri - 1) {-B+A^_ + 2B_A+A_ + BlB+) /2 
4/F^(r7 - 1) \b+A^_ - 2B_A+A_ - BlB+) f + 4:VP^{r] - 1) (A+Al + 2B_B+A_ - B'iA+) /2 
y/P^{r] - 1) {B+A^_ - 2B_A+A_ - B1B+) f + y/P^{r] - 1) {A+A'i + 2B_B+A_ - BlA+) /2 ) 

which is to be understood as a 4 x 3 matrix, and its corresponding analog. 
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